function resv = diffParamMatrix( ingConc, time )
    S = load('workParams.mat');

    sz = size(S.k);
    nsz = sz(2);
    
    res = zeros(nsz, 14);
    for i = 1:nsz
        res(i,:) = diffByParam(S.pf, S.pdm, cell2mat(S.k(i)), ingConc, time);
    end
    
    fprintf('\n');

    resv = res;
end

